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Abstract 

We analyze single transverse-spin asymmetries for hadronic pion production 
at large transverse momenta using QCD factorization. In the large xp region, 
leading contributions to the asymmetries are naturally produced by twist-3 
parton correlation functions that couple quark fields and gluon field strengths. 
With a simple model for these matrix elements, leading-order asymmetries cal- 
culated from QCD are consistent with data on pion production from Fermilab, 
and can be used to predict single-spin asymmetries at RHIC. We argue that 
our perturbative calculation for the asymmetries is relevant to pion transverse 
momenta as low as a few GeV. 
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I. INTRODUCTION 



Perturbative Quantum Chromodynamics (QCD) has been successful in interpreting and 
predicting spin-averaged scattering cross sections at large momentum transfer. Since quarks 
and gluons carry spin, we expect QCD to apply to hard spin-dependent scattering as well. 
However, high energy experiments with polarized beam and/or target have provided many 
theoretical challenges. For example, data on the spin asymmetries in deep-inelastic scatter- 
ing (DIS) of polarized leptons on polarized hadrons sparked a wave of theoretical effort 
in understanding the nature of the nucleon's spin [0. 

A spin asymmetry is the difference of two polarized cross sections, with opposite direc- 
tions of polarization, divided by their sum. Asymmetries can be obtained with both beams 
(or beam and target) polarized or only one beam (or target) polarized. The former is a 
double spin asymmetry, and the latter a single spin asymmetry. Depending on the direc- 
tion of the polarization, we can have longitudinal-spin asymmetries, if the polarization is 
along the beam direction, and/or transverse-spin asymmetries, when the spin is polarized 
perpendicular to the beam direction. 

Because of parity and time-reversal invariance, single longitudinal-spin asymmetries 
for single-particle inclusive production vanish for the strong interactions. However, 
experimentally-significant single transverse-spin asymmetries have been observed in A pro- 
duction, as well as pion production, for almost twenty years These single transverse- 
spin asymmetries are of the order of ten or more percent of the unpolarized cross section. 
Experimental results on pion production have been very consistent, and the effects persist to 
pion transverse momenta of several GeV, into the hard-scattering region, where perturbative 
QCD (pQCD) has had success in describing spin-averaged cross sections 0. The extension 
of the pQCD formalism to spin-dependent cross sections, however, has not been completely 
straightforward. It was pointed out long ago that QCD perturbation theory predicts 
vanishing single transverse-spin asymmetries at high px- Efremov and Teryaev later pointed 
out that a nonvanishing single transverse-spin asymmetry can be obtained in pQCD if one 
goes beyond the leading power 0J^. However, the relatively large size and peaking in the 
forward direction of observed effects remained a difficulty 0. 

Some time ago, using the example of hadronic direct photon production [|10| , we demon- 
strated that single transverse-spin asymmetries can be consistently evaluated in terms of 
generalized factorization theorems in perturbative QCD [|lT|. The asymmetries are presented 
as a sum of terms, each of which consists of a convolution of a twist-2 parton distribution 
from the unpolarized hadron, a twist-3 quark-gluon correlation function from the polarized 
hadron, and a short-distance partonic hard part calculable in perturbative QCD. The twist-3 
quark-gluon correlation functions reflect the interaction of quarks with the color field of the 



hadron |T0|,|T2|,|T3[ . In order to test this formalism, we need to have more than one physical 



process to extract information on these new and fundamental correlation functions, and to 



test their universality. Recent work has explored their role in the Drell-Yan process [jT4 
In this paper, we will not explore the physical interpretation of the correlation functions 
beyond what is currently in the literature. Rather, we concentrate on the extension of the 
formalism to pion production. 

In the forward region for pion production, where xp is large, we shall argue that leading 
contributions to the asymmetry depend on only one twist-3 matrix element (given in Eq. 
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(p5|) below), which couples two quark fields and one gluon field strength. This is the same 
matrix element that gives the leading contribution to single transverse-spin asymmetries in 



direct photon production |T0|. With a simple model for this twist-3 matrix element, we show 
that significant asymmetries can be generated, and that the asymmetries increase naturally 
as a function of a;^. Our simple model has two parameters: one for the normalization, 
and the other for the relative sign between the up and down quark correlation functions. 
Extrapolating from measured single transverse-spin asymmetries in vr"*" and tt" production 
in proton (I) -proton collisions [^], we fix these two parameters in our model. We can then 
derive both the sign and shape of the asymmetries for tt" production, as well as pion pro- 
duction in collisions with a polarized antiproton beam. Our results are consistent with data 
from Fermilab experiments. The model then predicts the normalization, xp and transverse- 
momentum dependence of the asymmetries at higher energies. These predictions can be 
tested at RHIC. 

The naive expectation for the twist-3 asymmetry, A^v, is A/^t, with A a nonperturbative 
scale from the twist-3 matrix element and ^T the transverse momentum of the observed 
particle. A pure 1/£t dependence, however, decreases quickly as £t increases, and becomes 
ill-defined when It is small. Consequently, one might worry that the range of It where the 
asymmetry is not too small, while ^T is large enough to use pQCD, is very limited, and 
that the region to study twist-3 physics might be too limited to be interesting. In fact, we 
shall see below that single transverse-spin asymmetries are a very good observable to study 
twist-3 physics perturbatively. 

In contrast to the naive expectation, for the kinematics of the Fermilab data, the \/^T 
contribution to is not the dominant source of the asymmetry. From dimensional anal- 
ysis alone, the asymmetry A^r admits two types of contributions, which are proportional 
to X£t/{—U) as well as XItH—T) ~ A/£t; with U and T Mandelstam variables. Their 
relative contributions can be determined by perturbative calculation. For large xp-, where 
the asymmetry is large experimentally, U is larger than T, but we shall show in this paper 
that the coefficient for the Mt/{—U) term is much larger than that of the A/£t term in this 
region (see Eq. ([76|) ). As we will see in Sec. 0, the transverse momentum dependence of the 
asymmetry is actually quite mild for It from less than 2 up to 6 GeV aX xp = 0.4, where 
much of the Fermilab data were collected. This conclusion is very encouraging for future 
applications of perturbative QCD beyond the level of leading twist. 

Our method and results can be generalized to single transverse-spin asymmetries in other 
single particle production. With the extracted information on twist-3 matrix elements, we 
can predict both the sign and magnitude of single transverse-spin asymmetries for any 
inclusive single-particle production, such as for direct photons, kaons, or other hadrons. 

Related work on single-spin asymmetries involves the incorporation of parton transverse 
momenta, either in parton distributions [^,0 or fragmentation functions ||T7HT9|. There is 
considerable evidence that at transverse momenta in the range of a few GeV, "fc^^-smearing" 
effects can be important in spin-averaged cross sections. It would seem natural to 

include them in the same range for single-spin asymmetries as well. Whether they should 
be thought of as the dynamical source of the asymmetry remains to be seen. The frag- 
mentation analysis requires the introduction of "chiral-odd" distribution functions |T7| , p!9 



which combine with the leading-twist transversity function [^,0 to produce nonvanishing 
asymmetries. We shall discuss how these effects can arise in the context of twist-3 factor- 
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ization theorems, but our explicit models will be based for simplicity on chiral-even parton 
distributions only. 

The twist-three analysis described here is in some sense a minimalist approach, depending 
on only light-cone variables, which we hope can serve as a benchmark for models which 
include models of both light-cone and transverse degrees of freedom. Other descriptions of 



single-spin asymmetries are based on multiquark interactions and orbital motion p4 
Interesting comparisons of different approaches may be found in and [P^]. 

Our paper is organized as follows. In Sec. |I|, we define single transverse-spin asymmetries 
in single particle production in hadronic collisions. We introduce generalized factorization 
formulas for the asymmetries, identify terms that we expect to dominate in the laige-xp 
region, review the factorization procedure at twist-three and leading order, and recall the 
leading-order spin-averaged cross sections to which we compare. In Sec. we present 
our explicit calculations of single transverse-spin asymmetries in hadronic pion production. 
We express these asymmetries in terms of short- distance partonic cross sections (coefficient 
functions), calculated in perturbative QCD, and non-perturbative twist-3 matrix elements. 
Using a simple model for the twist-3 matrix elements, we compare our calculated asymme- 
tries with experimental data in Sec. Finally, in Sec. 0, we give a summary of our results, 
and an outlook for the subject. We have also included an appendix, in which we review 
the application of parity and time reversal symmetry, and identify the list of chiral-even 
and chiral-odd twist-3 distributions and fragmentation functions that can contribute to the 
single-spin asymmetry for pion production. 



II. SINGLE TRANSVERSE-SPIN ASYMMETRIES 

A. Definition and General Considerations 

Single spin asymmetries are introduced for reactions in which only one particle is polar- 
ized. For example, consider single-particle inclusive production in a high energy collision, 

A{P,s) + BiP')^Cii)+X , (1) 

where A and B are the initial particles, with A polarized, where C is the observed particle 
(say, a pion) of momentum i, and where X represents all other particles in the final state. In 
order to fix the kinematics, we choose the center of mass frame of the incoming hadrons, with 
the 2;-axis along the momentum of the polarized hadron. We introduce two four-vectors, 
and n^, 

= (n+,n-,nT) = (1,0, Ot) , 
n'^ = (0,l,OT) , (2) 

with = = n"^, and n ■ n = 1. The incoming hadrons' momenta are P'* ~ \J and 
p'M ^ \/*Sy2, respectively. Invariants at the hadron level are defined as 

S = {P + P'f ^2P-P' 

T ={P-lf ^ -2P ■ i (3) 
U={P' - iY ^ -2P' ■ £, 
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where hadron masses are neglected. Given Eq. (j^), we next introduce 

21, T-U 



Vs s 



^T-^. (4) 

We now introduce a{i,s) as the cross section of the process given in Eq. ([l|). The 
spin-averaged cross section for single-particle inclusive production may be represented as 

a{i)^^[a{£,s)+a{i,-s)] , (5) 

and the corresponding spin-dependent cross section as 

Aa{i,s)^^[aii,s)-aii,-s)] . (6) 

The single spin asymmetry is often defined as a dimensionless ratio of spin-dependent and 
spin-averaged cross sections, 

A single longitudinal-spin asymmetry is denoted as Al, and a single transverse-spin asym- 
metry as An. We shall be concerned in this paper with An- For differential cross sections, 
the asymmetry can be defined as 

where Egd^a/d^i and E^d^Aa/d^i are the Lorentz invariant spin-averaged and spin- 
dependent cross section, respectively. In this paper, we will concentrate on single transverse- 
spin asymmetries in the forward region (i.e., large xp) where the asymmetries are largest 

i- 

Due to the symmetries of fundamental interactions, it is possible to have a vanishing 
single transverse-spin asymmetry, even though the corresponding total cross section a{i, s) 
itself is finite. For example, it was pointed out by Christ and Lee over 30 years ago ||2^ that 
time-reversal invariance forbids single transverse-spin asymmetries in inclusive deep- inelastic 
scattering (DIS) to lowest order in a^M- Let us review the reason. 

Consider a general inclusive lepton-hadron deep- inelastic scattering, which is the analog 
ofEq. d), 

L{i) + H{P, St) L{i') +X , (9) 

where L[i) and L{£') are unpolarized incoming and outgoing leptons of momenta, £ and £', 
respectively, and H{P, st) represents the polarized target hadron with its spin st perpen- 
dicular to the beam momentum. In the approximation of one-photon exchange, as shown in 
Fig. D, the inclusive cross section cr(sr) can be expressed as 
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aisT)(xL^''W^,isT) , (10) 

where the leptonic tensor, L^*^, is symmetric, and the hadronic tensor is given in terms of 
matrix elements of electromagnetic currents, 

W,,isT) oc {P,ST\jliO)Uy) \P,st) . (11) 

Applying parity and time-reversal (PT) and translation invariance to the matrix element in 
Eq. (|ll]), we obtain following relation. 



{P,st\jI{0)jM \P,st) = {P,-ST\jmjM \P^-^t) . (12) 

Combining Eqs. (pT]) and ([T^) , we find 

W^,{st) = W,,{-st) . (13) 

/^From Eq. (|^), we obtain the spin- dependent cross section for inclusive deep- inelastic scat- 
tering, 

Aa{sT) oc L'^'^ [W^u{st) - W^,{~st)] 

= , (14) 



where in the second line we use Eq. (|I3|) and in the third the symmetry of L^'^ when the 
lepton is unpolarized. From Eqs. (|^ and (0), it is clear that the single transverse-spin 
asymmetry for inclusive deep-inelastic scattering, A^^, vanishes to lowest order in oem- 

In hadron-hadron scattering, in contrast, the presence of multiple (initial-state or 
final-state) interactions prevents a simple decomposition like Eq. (|l3), and allows single 
transverse-spin asymmetries for final-state photons as well as hadrons [p|,|lO|. Experimen- 



tally, data from Fermilab show large single transverse-spin asymmetries in single pion pro- 
duction [Q, and at the same time, show no apparent single transverse-spin asymmetries in 
prompt photon production in the central (low Xp) region [ pSj] . 

Experiments at Fermilab for pion {71^,tt^) and prompt photon production were carried 
out with a 200 GeV polarized proton (or antiproton) beam on an unpolarized proton target. 
The observed single transverse-spin asymmetries of inclusive single pion production can be as 
large as 20 to 30% in the forward region. In addition to the large values of the asymmetries, 
a number of other interesting features are evident in the data. For example, a strong rise of 
the asymmetries with xp was observed for all pion charges. When the beam was switched 
from polarized proton to polarized antiproton, the same sign of the asymmetry was observed 
for 7r°, while the sign of the asymmetry for vr''", as well as vr", changed. Both beams had 
opposite signs of the asymmetries of 7r+ and vr^. 

Perturbative QCD was first used to study the effects of single transverse-spin asymme- 
tries by Kane, Pumplin, and Repko (KPR) 0. KPR calculated the single transverse-spin 
asymmetry for single hadron (pion) production in terms of a QCD parton model. By cal- 
culating the quark-quark scattering diagrams shown in Fig. |^, KPR found that the nonva- 
nishing single transverse-spin asymmetry for large-p^ reactions is proportional to the quark 
mass: Ajy oc ~ mg{P, sxlipTiplP, st) , where, for example, F = 7"*" 7577- Consequently, 
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the asymmetry vanishes in the scahng hmit (m^ — * 0). Although this calculation does not 
explain the observed large single transverse-spin asymmetries the fact that the result 
is proportional to the quark mass indicates that the single transverse-spin asymmetry is a 
twist-3 effect in QCD perturbation theory ||8| JT0| . |29[ . 



QCD dynamics, however, is much richer than the parton model. In addition to the par- 
ton mass effects just discussed, there are other twist-3 contributions. Because quarks are not 
exactly parallel to the incoming hadron beam, twist-3 contributions also arise from "intrin- 
sic" transverse momentum, which is proportional to T^^ ~ (P, StIiP'^^OtiPIP, St)- In addition, 
there are twist-3 contributions from the interference between a quark state and a quark-gluon 
state, which is proportional to T^j. ~ {P, st\'4'TA'j'iP\P, st) ■ Due to gauge invariance, T^j, 
and Taj. are not independent, and can be combined to form Tjj^ ~ {P,st\4'TDt4'\P,st), 
and/or Tp ~ {P,st\4'TF~^t4'\P,st), with F+t oc [D'^,Dt], where is the covariant 
derivative. Therefore, in addition to parton mass effects, single transverse-spin asymmetries 
can be proportional to the twist-3 matrix elements T^j. and Tp P, pil|P^ . These twist-3 
matrix elements involve three field operators {ipTDTip, ipTF^f'^^, or with the quark fields 
replaced by gluon field strengths pO]). Also, different choices for the Dirac matrices F in 



the operators give different twist-3 matrix elements (see the Appendix) |10 |. 

Because of their odd numbers of field operators, three-field twist-3 matrix elements do 
not have the probability interpretation of parton distributions, which are proportional to 
matrix elements of twist-2 operators, ipTil) or F^j-F^t- In principle, however, they are as 
fundamental as the parton distributions. Measurements of twist-3 distributions, or three- 
field correlation functions, provide us new opportunities to study QCD dynamics. 



B. Factorization and the Valence Quark Approximation 

As we have seen, spin-dependent asymmetries for hadronic pion production with one 
hadron transversely polarized vanish at large momentum transfer Nonvanishing values 
of the single transverse-spin asymmetry signal non-leading power contributions. According 
to the basic factorization theorems the leading power spin-averaged cross section for 



the production of a pion with large transverse momentum ip can be factorized into four 
separated functions, as sketched in Fig. |, 

(JA+B^tt = (Pa/A{x) (g) (pb/Bix') (g) CTa+b^c ® Dc^7t{z) , (15) 
abc 



where Y^abc represents the sum over parton flavors: quark, antiquark and gluon. In Eq. (^5|), 
4>a/A{x) and (pb/six') are probability densities to find parton a of momentum xP in hadron 
A and parton b of momentum x'P' in hadron B, respectively. As noted above, they may be 
interpreted in terms of expectation values in the hadronic state of two-field matrix elements, 
for example ipTip or F^tF^t- Dc^t^{z) is the fragmentation function for a parton c of 
momentum pc = C-/ z to fragment into a pion of momentum £, and (Ja+b^c is a short- distance 
partonic part (the Born cross section plus corrections), calculable perturbatively order- by- 
order in Us- The symbol ® in Eq. (|T5|) represents the convolution over the corresponding 
parton momentum fraction. In terms of the Lorentz invariant differential cross section, 
Eq. (ITSD can be written as pi 
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Ee -pj-^ ^^<^«m(^) J dx (j)b/B{x') I — [ Ec — ) — • (16) 

abc 




The predictive power of Eq. ( p!6|) depends on independent measurements of the non- 
perturbative functions, (pa/A^ (pb/s and -Dc-^tt? and the calculation of the partonic part 
E^d^aa+b^c/d^Pc- 

Just as for most other physical observables calculated in perturbative QCD, the predictive 
power of the theory for twist three relies on factorization theorems [^. Physical observables 
that depend on the transverse polarization of a single hadron are typically power corrections 
to the total cross section, in comparison with spin-averaged or longitudinally polarized cross 



sections. In Ref. |^2[, for a physical observable with a large momentum transfer Q, we 
extended the factorization program to 0(1/Q^) corrections for unpolarized hadron-hadron 



cross sections, and in |[TT| , |33[| to 0{1/Q) corrections in polarized cross sections. 

Following the generalized factorization theorem |^T],^, the transverse spin-dependent 
cross section for large £t pions, Aa{sT), can be written in much the same way as the spin- 
averaged cross section, Eq. (|l^), as a sum of three generic higher- twist contributions, each 
of which can also be factorized into four functions, 

AaA+B^nisr) = ^i/A^^l' ^2, St) <^ (f)b/B{x') ® Ha+b^c{sT) ® Dc^„{z) 

abc 

+ E ^i/Al^, ^t) ® 0f/k^'l> 4) ® K+b^cisT) ® D^^^iz) (17) 



abc 



abc 

+ higher power corrections , 

where Y^abc represents sums over parton flavors: quark, antiquark and gluon, and where 
4>b/B{x') and Dc-,tt{z) are standard twist-two parton distributions and fragmentation func- 
tions, respectively. In Eq. (|17|), the first term corresponds to the process sketched in Fig. §a, 
and the second and third terms correspond to the ones sketched in Fig. 

For the first term in Eq. ([T7|) , nonvanishing contributions to A(t(st) come from twist-3 



parton distributions (correlation functions) 0a/^^(a;i, X2, st) in the polarized hadron. For the 
second and third terms, the contributions to A(j(st) involve the twist-2 transversity distri- 



butions 6q^\{x, St) pT| , p2| . Because the operator in the transversity distribution requires 



an even number of 7-matrices |^ , ^ , the second term and third terms in Eq. ([TtD also 

fr>\ 

include a twist-3, chiral-odd parton distribution 06/5 (a^'i, 3^2) from the unpolarized hadron 
B, or a twist-3, chiral-odd fragmentation function, D^J^^{zi, Z2). In the factorized form of 
Eq. ([l^), PT invariance may be applied in a manner analogous to the treatment of the DIS 
cross section given above. In this case, however, PT invariance allows nonzero An for a 
limited number of functions, as discussed in the Appendix. 

As in the spin-averaged cross section, Eq. (p!5D, the hard-scattering functions Ha+b^c{sT) 
are the only factors in Eq. (0) that are calculable in QCD perturbation theory. The 
calculation of the if's depends on the explicit definitions of the twist-3 distributions, for 
example 0^/^^(xi, 0:2, sr), and the predictive power of Eq. ([T7|) rehes on the universality of 
the new twist-3 distributions [|Tl|,B31. 
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Eq. ([ITD illustrates the typical complexity of higher-twist analysis: even at first nonlead- 
ing twist, whole new classes of functions begin to contribute. This complexity is particularly 
difficult to sort out for physical observables to which leading-twist terms contribute. The 
combination of small effects and complex parameterizations has made the extraction of 
higher twist distributions from the data difficult, despite the considerable effort that has 
been invested in the formalism 0,^. Q The vanishing of single-spin asymmetries at lead- 
ing power solves one of these problems, the masking of higher twist by leading twist. Beyond 
this, however, it is clear that to fully disentangle all of the functions contributing to Eq. (|1^ 
would require a constellation of data and a level of analysis far beyond what is currently 
available. Turning specifically to the first term in Eq. (p^Tf), we observe that the index a refers 
to pairs of partons, and that the functions 0^^^ are correspondingly functions of two mo- 
mentum fractions. In addition, even assuming that we knew this set of functions, we would 
still be faced with the chiral-odd distributions and fragmentation functions in the second 
and third sums in Eq. (p^Tf ). We would like to suggest, however, that by restricting ourselves 
to the limited kinematic range of large xp for the observed particle, we may simplify the 
analysis greatly, and construct a simple model that explains the available data, and that 
provides extrapolations to higher energies and momentum transfers. 

We are going to present a calculation of the Xaxge-xp asymmetry at moderate or large 
^T■, in terms of the chiral even functions ^^^^^(xi, 0:2, St) only (first line of Eq. (|T^). In these 
functions, we will consider only combinations of valence quark fiavors with gluons. We will 
not find it necessary to specify these functions for all values of Xi and X2, but only for the 
line Xi = X2i at which the gluon carries vanishingly small momentum fraction. We will 
refer to this set of simplifications as the valence quark-soft gluon approximation below. In 
this model, we thus neglect potential contributions from the transversity, coupled with the 
chiral-odd twist-3 distributions and fragmentation functions identified in the Appendix. We 
hope to explore these contributions elsewhere, but in the absence of independent information 
on the transversity, it seems natural to test the plausibility of a model based on chiral-even 
distributions alone. 

First, consider our restriction to valence quarks. Given that single transverse-spin asym- 
metries were measured at Fermilab with a 200 GeV polarized beam |Q, only partons (a and 
h in Eq. (0)) with large momentum fractions will be relevant for large xp or ^J•■ Because 
parton-to-pion fragmentation functions vanish as z 1, the effective momentum of the 
fragmenting parton, pc = i/z, should be much larger than the pion momentum i. There- 
fore, the dominant contribution to the cross sections in the central region should come from 
X ^ x' in Eq. (|13), with, in addition, x much larger than xt ~ 0.25, which corresponds to 
£t ^ 4 GeV at -Ebeam = 200 GeV. In our calculation we will concentrate on the forward 
region, where xp is large. Similarly, in this region the dominant contributions to the cross 
section come from x considerably larger than xt (i.e., x > 0.25) even for relatively small ip- 
For large x, there are few gluons or sea quarks from the beam hadron. Therefore, in our 
numerical calculations, we will keep only valence quarks from the polarized beam. That is, 
J2a in Eq. (|17D now runs over only up and down valence quarks, coupled with a single gluon 



^We may note recent progress based on models of higher twist in deeply inelastic scattering and 
fragmentation inspired by renormalon analysis |3f:-38|. 



field. 

In presenting this argument, we are well aware that in principle the flavor content of the 
twist three distributions may be totally different than those of twist two. Nevertheless, we 
consider it by far more natural to assume that three-field correlations at large Xi will be 
dominated by the same flavors as in the two-field, parton distribution, case. We recognize 
that this remains, however, an assumption. In our case, it means that we shall keep only 
valence quarks from the polarized beam, accompanied in twist-3 by gluons. In particular, 
we shall not consider three-gluon matrix elements . 



We now turn to the question of "soft gluons" . To anticipate, the twist-three asymmetry 
involves only two classes of contributions in Habc{xi,X2)- One of these is proportional to 
6{xi — X2), and the other to S{xi), i = 1,2. The first case sets the momentum carried 
by the gluon field in the twist-3 matrix element into the hard scattering to zero, leaving 
the momenta carried by the two quark fields in the combination ipF'^'^ip diagonal^. In the 
second, one of the quark fields {ip or ip) carries vanishing momenta. We refer to these two 
possibihties as "soft gluon" and "soft fermions" poles, respectively ||10| . Soft gluon terms are 



(3) 

typically accompanied by derivatives of the parton functions (/'^/^, while soft fermion terms 



are not. We have emphasized in Ref. |T0[ that terms that involve derivatives with respect to 
distributions tend to be strongly enhanced near the edges of phase space, relative to those 
without derivatives. We shall see this in our explicit model below. We shall assume, in 
fact, that it is this effect that is primarily responsible for the experimentally-observed rise 
in single-spin asymmetries toward xp = 1- We therefore suggest that only terms in which 
such derivatives occur need be kept, in order to describe the laige-xp single-spin asymmetry. 
In summary, only soft-gluon terms, from the first line of Eq. ([T7|) produce the shape of the 
large asymmetries observed in the data in the forward region, and for these terms xi = X2- 
To set the stage for the explicit calculations of the next section, we first give an example of 
leading-order factorization at twist three for the spin-dependent cross section, following the 
method of Ref. |jlO| . This will enable us to trace the origin of twist-three spin distributions. 



and of the poles that underly the valence quark-soft gluon approximation that we have just 
described. 



C. Twist-3 Factorization at Leading Order 

The twist-3 correlation functions, X2, Sr), depend on two parton momentum 

fractions, while twist-2 parton distributions, which are probability densities, depend on 
only one. Considering the effort and data needed to determine the parton distributions, 
it appears a difficult task to get a full description of these twist-3 distributions. From the 
general Feynman graphs contributing to the if's in Eq. ([T^), as shown in Fig. ^ it is also 
clear that there are many diagrams, even at lowest order. Their treatment is simplified, 
however, by taking advantage of the relation of the asymmetry to the pole structure of H 
p!0| . This will enable us to evaluate Acr(sr) in Eq. (|T^ efficiently. Indeed, we will find that 
Aa^sx) depends on the twist-3 distributions through only a single independent momentum 



^Note, there are no "soft gluons" in the short-distance functions H. 
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fraction, with the other fraction fixed by a pole. To see how this comes about, we consider 
a specific set of contributions, associated with the three classes of diagrams shown in Fig. p. 



We first discuss the analysis of these diagrams according to the method of Ref. |T0[, and 
then briefly discuss other possibilities, reviewing why we expect those of Fig. ^ to dominate 
the asymmetry in the large-x^ region. 

In our valence quark-soft gluon approximation, introduced in the last subsection, the 
fermion flavor a from the polarized hadron in Figs. ^ and |^ runs over valence quarks only, 
while parton h from the unpolarized hadron can be a gluon, valence quark or sea quark. 
We start from these three classes of diagrams, and derive below the factorized form for 
the spin-dependent cross section Acr(sV'). The hard-scattering diagrams of Fig. |^ are all 
embedded in the overall process shown in Fig. 0a. The top part of this general diagram 
is proportional to the expectation value of an operator of the form ipA^ip in the polarized 
incoming hadron state \P, st), while the bottom part includes the hard subprocess, as well 
as the target hadron matrix element and the final-state pion fragmentation function. In 
Fig. ^ ki and k2 are valence quark momenta, and a is the Lorentz index for the gluon field. 



We work, as in Ref. [|T^, in Feynman gauge. To derive a factorized expression for these 
contributions, we must separate spinor and color traces, as well as sums over vector Lorentz 
indices between the functions T and S. 

After separation of all traces by a Fierz projection (see Appendix), the two functions T 
and 5* are connected only by the two momentum integrals that they share. The leading con- 
tributions of the general diagram shown in Fig. ^ can then be represented by the factorized 
diagrams shown in Fig. |^, and can be written as 

dAaisT) = ^ E / ^ ^ lUh, h, st) S^{k„ k,)] , (18) 

where 1/25* is a flux factor, J2a I'uns over only valence flavors, Ta{ki, ^2, st) is proportional 
to the matrix element of the operator, (27r)[?/'a7"'"y4"'"'?/'a]/2P"'"^, and Sa{ki, ^2) represents the 
bottom part of the general diagram shown in Fig. |^. The function Sa{ki, ^2) is contracted 
with [(1/2)7 ■ -P-Pcr] Ca/(2vr), where the factor (27r) is due to the normalization of twist- 
3 matrix element T, which we will specify below. The color factor Ca is left from the 



factorization of color traces between Ta{ki,k2,ST) and Sa{ki,k2) ||T0[- With the function 



Ta{ki^ k2, St) oc i/jAi/j^ the corresponding Ca is defined for all valence flavors a as 

(C^f = {j^) it%j , (19) 

with = 3 colors, B the gluon color index, and with quark color indices ij. The matrix 
{t^)ij is the SU(3) generator in the defining representation of the group. 

The next step in the factorization procedure is the "coUinear" expansion which 



will enable us to reduce the four-dimensional integrals in Eq. (0) to convolutions in the 
momentum fractions of partons, as in Eq. (|T7|). Expanding Sa in the partonic momenta, ki 
and ^2, around ki = xiP and k2 = X2P, respectively, we have 



Saih, h) = Sa{xi, X2) + -^{xi, X2) {ki - XiPy 



dkl 

dk^o 



+ ^ix^,X2)ik2-X2P)'+... . (20) 



11 



This expansion, substituted in Eq. ([T8|), allows us to integrate over three of the four com- 
ponents of each of the loop momenta ki. The top part of the diagram then becomes a 
twist-three light cone matrix element, convoluted with the terms of Eq. (|20| ) in the remaining 
fractional momentum variables x,-. 



As stressed in Refs. P,p^, some of the matrix elements that result from the collinear 



expansion can have nontrivial spin-dependence. It is at this stage that the pole structure 



of the hard scattering begins to play an important role. In fact, as shown in Refs. PJTI] 
and below, nonzero spin dependence is found only from pole terms in the hard scattering. 
Without these poles, the symmetries of the strong interaction force the asymmetry to vanish, 
in much the same fashion as for DIS above. Indeed, the poles provide exactly the sort of 
multiple interactions that are absent in DIS at lowest order in QED. The first term in 
the expansion, Eq. (|20|) , Sa{xi,X2), does not contribute to A(t(st) when combined with 
Ta{ki,k2,ST) in Eq. ([ISD , because it lacks true initial- or final-state interactions. We will 
therefore drop it below. 

Let us next look for poles in the diagrams of Fig. ^ from the remaining terms in Eq. ( pOD , 
and identify the relevant twist-three matrix element. All of the diagrams in Fig. ^ provide 
a pole at xi = X2 when ki = XiP{i = 1,2). As we will show below, these poles have the 
property that 

^{XuX2) = -^{XuX2), (21) 

for Xi = X2- This equality is to be interpreted in the sense of distributions, since Sa is singular 
at Xi = X2- Substituting Eq. (^) into Eq. ( pOD and neglecting higher order derivatives, we 
have 



Saiki,k2) ^ T^{x,,X2) [u;f^{k2 - k^),] , (22) 



dSa 

dkl^ 

where the projection operator u^"' is defined as u^"' = g^"' — -n^n^. Substituting Eq. (^2]) 
into Eq. (|18|) and performing the integration over the non-longitudinal components of the 



fc's, we derive 



dAaisr) = J dxidx2 



T^^J{x^,X2), (23) 



where the integration over xi (or X2) will be fixed by the corresponding pole in dSa/dk2, 
and where e^'^Tnn jg fjgg^ed as 

^psTun ^ ^papv _ ^24) 

The function Tp^\xi, X2) for flavor a in Eq. ( P^D is one of the twist-3 distributions introduced 
in Ref. 0, 



T?'J{xuX2) = J ^f^LMe-i^+?/r+K-2-xi)P+y,- 



x{P,ST\MOh+ e™F/(y2-) Myi)\P,ST) . (25) 
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The ordered exponentials of the gauge field that make this matrix element gauge invariant 
have been suppressed |32,33|. It is easy to show that is real. Parity ensures that 
^pnns^ and time reversal invariance then implies that it is an even function of xi and 



"-Fa 

X2, 



Tj^\xi,X2) =T^y{x2,xi) . 



iV), 



(26) 



These properties are valuable in isolating nonvanishing asymmetries. For instance, the fact 
that Tp^^ is real ensures that only the poles of S in Eq. ( ^5] ) can contribute. 

Having factorized the twist-3 distribution we now factorize the remaining func- 

tion [ie'"^^"'"'dSa/dk2] in Eq. (|2^) into a perturbatively calculable partonic part Ha+b^c, a 

. At the 
and can 



corresponding target parton distribution (p^/B and a fragmentation function De- 
leading power, diagrams contributing to Sa{ki, /C2) can be represented as in Fig. | 
be factorized as 

'^''''■[Ma+biki,k2,k')B,ik',P')] 



Sa{k 



1,^2) 



E 

b 

E 



(2vr)4 
dx' 



X' 



Ma+b{ki,k2,x')(l)b/B{x') 



(27) 



where J2b runs over all parton flavors, and 4>b/B{x') is a twist-2 parton distribution for flavor b, 
for the unpolarized target hadron B. We use the matrix element definitions of twist-2 parton 
distributions given in Ref. [^. Similarly, as shown in Fig. ||b, the factor Ma+b{ki,k2,x') 
in Eq. (p7|) can be further factorized into a convolution of Feynman diagrams, calculable in 



(28) 



perturbation theory, with standard twist-2 fragmentation functions, 

Ma+biki, k2, x') ^ J dz Ha+b^ciki,k2,x',pc = i/z) Dc-*^{z = i/pc 



where the Ha+b^c are given by the diagrams of Fig. ^ The fragmentation functions Dc->-k{z = 
i/pc) are also defined as matrix elements in Ref. i39(| . 

Finally, substituting Eq. (|2|) into Eq. (0), and Eq. into Eq. (0), we derive a 
factorized expression for Aa{sT) in the form of Eq. (ITT 



dAa{sT) = ^11 j dzD^^^iz) J 

abc 



dx' 
x' 



dxidx2T^\xi, X2) 



X 



le 



d_ 



Ha+b^c{ki = XiP, k2 = X2P, x',pc 



(29) 



where the integration over either xi or X2 can be done by using the pole in Ha+b^c- This 
results in a factorization with only a single momentum fraction for each of the incoming 
hadrons, similar to that for the spin-averaged cross section in Eq. (0), with (pa/A^x) replaced 

by Tj^\x,x). In order to use this factorized formula for single transverse-spin asymmetries 
in pion production, in the following section we will evaluate the diagrams shown in Fig. ^ 
with off-shell momenta ki and k2- In each case, we will verify Eq. (0), or equivalently, 
observe that 



dH 



{Xi,X2 = Xi) 



dH 



(Xi,X2 = Xi) 



dk2, dki^ 
where, again, the equality is to be interpreted in terms of distributions. 



(30) 
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D. Leading Contributions in the Forward Region 



Before entering into the detailed calculations of the hard-scattering functions Ha+b^c in 
Eq. p9|), we return to issue of why we believe that the dominant contribution is given by 
the tP in Eq. (pQ]). We have already indicated that this is due to the derivative structure 
of these contributions. Let us see how these derivatives arise. 

/,From the diagrams shown in Fig. with the momenta Pc and x'P' fixed, we get four 
typical sources of ki {i = 1,2) dependence: (1) /cj-dependence in (5(L(/cj)^) with L the mo- 
mentum of the unobserved final-state parton, (2) /cj-dependence in the propagators which go 
on-shell when fci = /c2, (3) fcj-dependence in the off-shell propagators, and (4) fcj-dependence 
in the numerators. The derivatives of Ha-^b^c{ki,k2,x',pc) with respective to ki have the 
following features: 

1. {d / dki)5{L{kiY) gives 5'(L(xjP)^), and its contribution to /S.a{sT) is proportional to 
{d/dx)Tp'\x, x) after integration by parts; 

2. (d/dki) on a propagator that is potentially on-shell changes a single pole to a double 
pole, and the resulting integration over the double pole makes the contribution to 
Acr(sr) proportional to {d/dx)Tp'\x,x); 

3. (d/dki) on an off-shell propagator does not change the pole structure, and its contri- 
bution to Aa{sT) is proportional to tP{x,x) without a derivative; 

4. (d/dki) on fcj-dependence in the numerator gives contributions to A(x{st) proportional 
to tP (x, x) without derivatives. 

As we have pointed out earlier, we are interested in the asymmetries in the forward region, 
where xp is large. Asymmetries in this region are dominated by large net momentum fraction 
X from the polarized beam parton, coupled with relatively small momentum fraction x' from 
the partons of the unpolarized target hadron. Since all distributions vanish as a power for 
large x, as (1 — x)^ with l3 > 0, {d/dx)TP{x,x) ^ tP{x,x) when x — > 1. Therefore, 
in the forward region, terms proportional to derivative of the distributions dominate. 
In order to simplify our calculations of the largest effect, we keep only these terms. Thus, 
in Sec. |T| we will keep only those contributions corresponding to items (1) and (2) listed 
above. 

Turning, finally, to other possible contributions in Eq. (|1^, we observe that it is only 



the matrix element that inherits derivative terms, as a result of the collinear expansion 



involving soft gluon poles. Soft fermion poles, of the sort discussed in Refs. p|,|TO| have no 



such derivatives at leading order. Soft-fermion poles also do not correspond to the valence 
quark approximation identified above, since they require one of the quark fields to carry zero 
momentum fraction. These features of the calculation follow exactly the same pattern as 



for direct photon production, as treated in Ref. ||T0[, and we shall not repeat them here. It 
is only necessary to emphasize that the contributions from Fig. |^ are the complete set 
of derivative contributions at twist three and leading order, for the first (chiral even) term 
in Eq. 
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E. Spin- Averaged Cross Sections for Hadronic Pion Production 



In order to evaluate the asymmetries, defined in Eq. (|), we need to compute the leading- 
order spin-averaged cross section. QCD perturbation theory has been generally successful 
with experimental data on spin-averaged cross sections for inclusive single-pion production 
at large transverse momentum ||^. At leading order in a^, only 2^2 Feynman diagrams, 
shown in Fig. ^ contribute to Ecd^aa+b^c/d^Pc- In terms of scattering amplitudes, the 
leading order Ecd^aa+b-*c/ d^Pc can be expressed as 

^ ^ '5{s + i + u), (31) 



d^Pc 



IGvr^s 



M. 



where M is the spin-averaged amplitude. In Eq. (|31|), invariants at the parton level are 
given by 



s = {xP + x'P'f = xx'S , 
i = {xP — Pc)"^ = xT / z , 
u = {^x'P' -pc? = x'U /z , 



(32) 



where S, T and U are defined in Eq. (^ . 

In the valence quark approximation, using the 5-function in Eq. (31) to fix the x'- 
integration in Eq. (|T^), we find the spin-averaged cross section for pion production at leading 
order in a,, 



E, 



1 dx 



a _ a; ^ dz ^ 

d ^ iS CL C ■^min ^ "^niin OCiS \ / ^ 



dx' / , —xT/z 
'V y ~ xS + U/z^ 



X qa{x) 



G{x') CTag^c +^q{x')a, 



aq^c 



(33) 



where J2a runs over up and down valence quarks, and J2q over quarks and antiquarks. In 
Eq. (P^D, the integration limits 2min and Xmin, and variable x are given by 

-{T + U) 



S 
-U/z 



X 



S + T/z 
—X T / z 



(34) 



xS + U/z' 

where 5, T, U are defined in Eq. (|^), and xp and xp in Eq. (^). The short-distance partonic 
parts, CTag^c and daq^c, in Eq. (|3|), are given by ||] 

4 



(7, 



ag^c 



/ SU \ 4 / —U S \ S U 



(35a) 



cr, 



aq^c 



, 4 fs^ + u^\ ^ . 

Oar - I I + Oqc- 



4 (s^ + F 
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t2 



9 



4 fP + 



(35b) 
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where s,i,u are defined in Eq. (|32D . 

Since we are interested in the large xp region, we have i'^ ^ £^ and T <^ U < S. 
Therefore, leading contributions to the cross section given in Eq. (^) come from the t- 
channel diagrams (the first diagrams in Fig. |^a and Fig. ^), or equivalently, the 1/r term 
(i.e., first term) in Eq. (^5]a) and Eq. (|35| b). Consequently, for leading contributions in the 
forward region, incoming parton a has the same flavor as fragmenting parton c. Therefore, 
in the valence quark approximation, we keep only 0^^^^+ for vr"*" production; D^-^.^- for tt" 
production, although we keep both D^^^o and -D^^^o for 7r° production. 



III. CALCULATION OF THE ASYMMETRY 

In this section, we present our calculation of the single transverse-spin asymmetries in 
pion production in the valence quark-soft gluon approximation described in the previous 
section. We derive analytic expressions for the spin-dependent cross section, A(t(£, sr), 
which is needed to evaluate the asymmetries. 



A. Quark-Gluon Subprocesses with Initial-State Interactions 

Consider the two diagrams with poles from initial-state interactions, as shown in Fig. ^a. 
We parameterize the parton momenta ki as 

ki = xiP + /cij,, and k2 = X2P + /c2t , (36) 

with the fcjj, two-dimensional transverse momenta. The remaining momentum components 
do not enter at twist three. The pole in the diagram at the left of Fig. ^ is given in these 
terms by 

(37) 



{x'P' + k2 — ki)"^ + it {x2 — Xijx'S + {k2j, — ki^Y + ie 
The derivative of this pole with respect to k2j, (or ki^) vanishes as ki^ -^0. The diagram 



on the right has the same feature. Therefore, following the arguments of Sec. |IID| above, 
the leading contribution to A(t(st) in the diagrams in Fig. ^ is from the derivative of the 
phase space ^-function only. 

Let Li and L2 be the momenta of the unobserved partons in the diagrams to the left 
and right, respectively in Fig. |]a. We have 

Li = x'P' + XiP + ki^ - pc, L2 = x'P' + X2P + k2j, - Pc ■ (38) 

Taking the derivative with respect to k^ and k2, we obtain 

^^6{Ll) = {-2pJ6'{Ll) , (39a) 

^ 6{L\) = , (39b) 



dk 



p 
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d 
d 



5{Ll) 



5{Ll) 



0, 



(39c) 
(39d) 



In deriving these relations, we have used that p is a transverse index. After taking the 
derivative with respect to the ki on the 5-functions, we can set kj^rj, to zero in the remainder 
of each diagram. For the diagrams in Fig. the poles giving the leading contributions are 
from 



9s{-i) gpf3 



{x'P'),Pp - {x'P')pP,] 



[x'P' +{x2-xi) P]'^ + ie 



-1 



R{xi,X2) = Qs {-i)gp(3 



Xi — X2 — i€ 
-1 

X2 — Xi + ie 



(40a) 
(40b) 



where L and R represent the diagrams at left and at right, respectively. In Eq. (^OD, 
Qs = v^47ra7 is the strong coupling. In the following discussion, we absorb the overall {—i) 
in Eq. ( |40| ) into the color factor for the subprocess. Using the distribution identity 



1 



X2 — Xi ± ie 



P 



X2 - Xi 



=F in5{x2 - Xi] 



(41) 



for the poles in Eq. (|^) and keeping the imaginary contribution of the pole, we can express 
the contributions of the diagrams in Fig. ^ as 



d 

-^{Ha^{Xi,X2,x',Pc) + Ha^{XuX2,x',Pc)) 



9s 



2ttxi 



-H. 



2^2 



'x2,x',pc) i7T6{Xi - X2){2pcp) 5' {L^) 



d 

■-^{Ha^{Xi,X2,x',Pc) + Haj,{Xi,X2,x',Pc)) 



2nx 



-H2^2{xi,x',pc) iTT6{x2 - xi){2pc„) S'{L-^^) 



(42a) 



(42b) 



where subscripts and represent the left and right diagrams of Fig. ^ja. In Eq. (^2|), 
H2^2{xi,x',Pc) with i = 1,2 is proportional to the imaginary part of the 2^2 partonic 
forward scattering amplitude shown in Fig. llO, 



H2^2{Xi,x',Pc) 



167r2 



M 



a+g^c 



(43) 



where the matrix element squared, M^^^_^^ , is the same as that in Eq. (0), except for 
the color factor, C^, due to the extra initial-state interaction. Combining Eqs. (19) and 
(PH), the factor Cg is given by the color structure of the partonic diagrams shown in Fig. |^a, 
contracted with a common factor [(— ■i)2/(A^^ — 1)] {t^)ij, where B and ij are color indices for 
the gluon and quarks from the polarized hadron. The factor l/27r in Eq. ( ^2]) was explained 
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in the text following Eq. (|T8]), and the factor is due to the definition of H2^2{xi,x',pc), 
where incoming quark lines are contracted with (1/2)7- i^iP)- Eq. (^21) shows that Eq. (|30D 



is satisfied when ki,^ = 0. 



Substituting Eq. (0a) into Eq. (|29|) , we have a complete factorized form for the spin- 
dependent cross section from the diagrams shown in Fig. lOa, 



dH 



J ^D^^A^) J dx'G{x') J dxTP{x,x) U 



dAai^g^^isT) 



d^Pc 



(44) 



where the factor l/z"^ is due to the phase space difference between d^i/{27c)^2E£ and 
d^Pc/ {2'7iY2Ec, and the partonic hard part, EcdAal^_^g^^/ d^pc, is given by 



E, 



dAai^g_^^{sT) 



d'^Pc 



IGvr^s 



M 



a+g^c 



' 6'{s + i + u) 



(45) 



In Eqs. ( ^4]) and (|45|) , superscript / indicates the contribution from a partonic subprocess 
with an initial-state pole, and subscript g represents the quark-gluon subprocess. In deriving 
Eq. (p^ , we renamed the integration variable Xi in Eq. (|42|) as x. The factorized spin- 
dependent cross section given in Eq. ( |1^ ) is very similar to the factorized form for the spin- 
averaged cross section in Eq. (p!6|), with the unpolarized parton distribution 4>a/A{x) replaced 
by the twist-three correlation function Tp'\x, x). The partonic hard part in Eq. (^) is also 
very similar to that in Eq. (^). For the spin-dependent case, the derivative of the (5-function 
is just the derivative with respect to the parton momentum ki in Eq. (pO[), which comes 
from the collinear expansion. The factor e^TPcnn (|45|) is necessary for a nonvanishing 

asymmetry. 

After partial integration over x, we can reexpress the derivative of the ^-function as 



/ dx 5' is + i + u)F(x) = [ — 
J J x'l 



dx 



-5{s + i 



u] 



'S + T/z 

for any smooth function F{x). Using Eq. (|45|), we thus rewrite EidAa^lsT) / d^i as 



(46) 



E. 



dH 



1 dx 



St inn 



1 



x'S^Tjz^ 



in X xS + U / z 

d ( T^/^\x,x) J 



dj^ (,_^xriz^ 

x' \ xS + U/z 



-X 



dx 



X 



(47) 



where Zmm and Xmin are given in Eq. (|3J) , and S, T and U are defined in Eq. (H) . In Eq. 
the spin-dependent cross section E^dAagi^sx) I d^C has almost the same factorized form as 
the spin-averaged cross section shown in Eq. ([331) . The extra factor of 1/z is due to the 
replacement of p^ by C in the e-tensor of Eq. (^5)) . The dimension of l/{x'S + T/ z) due to 
the derivative of the 5-function is balanced by the dimension of ^ in the e-tensor and the 
dimension of the twist-three correlation function t'^\x, x). In our definition, the twist-three 

[s, £, m), in 



correlation function has the dimensions of energy. The partonic hard part, iff, 



ag- 



Eq. (^ plays the role of aag^c in Eq. (^). It is given by Cg |M„+g_,j2 in Eq. (^, which 



represents the 2 — > 2 matrix element squared in Eq. (|43|), but with a different color factor 
due to the extra initial-state interaction. 
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B. Quark-Gluon Subprocesses with Final-State Interactions 



The diagrams shown in Fig. ^ represent final-state interactions of the fragmenting par- 
ton. As with the contributions from initial-state interactions, these diagrams also have a 
derivative with respect to fcf and kl^ of the phase space 5-function associated with the un- 
observed final-state parton, of momentum Li or L2. Similarly to Eq. (^0|), the final-state 
poles giving leading contributions are given by, as sketched in Fig. |ll]. 



^9s{l- Pc 

R{xi,X2) ^ gs {1 -Pc 



7 ■ P7 ■ (pc + (a^i - X2)P) 
{Pc + {xi - X2)PY + ie 
1 



Xi 



X2 + ie 
1 



X2~ xi — le 



(48a) 
(48b) 



where the factor (7 ■ p^) "will be absorbed into the 2 — >■ 2 hard-scattering function. Similarly 
to Eq. (^71) , we obtain the contribution from the derivative of the 5-function for a final-state 
interaction. 



S 



E 



1 dz 



dx 



1 



X xS + U / z 



dx' f , —xT/z 
~ \ ~ xS + U/z^ 



STinn 



x'S + T/z^ 



Gix') 



9 /TfW) 
dx \ X 



-X- 



H. 



ag^c 



(49) 



where superscript F denotes the final-state interaction. The only difference between 
EidAa^{sT)/d^i in Eq. (H) and Eid/^alM / d^(. in Eq. (|^ is the color factors in the 
partonic hard parts. The hard part H^g_Js,i,u) in Eq. (^91) i-^ r<p'\]\j' 12 



in Eq. (H) is given by \M 



I 

a+g- 



which has the same kinematic dependence as H^g^^{s, t, u) in Eq. (^T]) , but a different color 
factor, Cf^, due to different color structures in final-state compared to initial-state inter- 



9 ' 



actions. Similarly to C^, is computed by contracting the matrix [2/(A^^ — 1)] {t^ 



Eq. (|TUp, into the diagrams. 

In addition to the contribution from the derivative of the (5-function, the diagrams shown 
in Fig. |b also give leading contributions, proportional to {d/dx)Tp'\x,x), from the double 
pole which results when the derivative (d/dki) acts on a propagator that goes on-shell at 
Xi = X2- Consider the final-state interaction in the diagram at the left in Fig. p!T|a. The pole 
giving the leading contribution is from the factor 



L{kir^, k2j 



9sl-Pc 
9sl-Pc 



7 ■ P7 • (pc + fci - ^2) 

{Pc + ki- ^2)2 + ie 

1 



xi- X2 + xo{ki.^, A;2j,) + ie 

7 • (feiT - hr)!- P 

2P ■ Pc [xi - X2 + Xo{ki^, k2j,) + ie] 



(50) 



where Xq is defined as 
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2{h 



Pc + {ki 



2P-Pc 
as ki and k2„ — > . 



(51) 



In deriving Eq. (|50[), we used the parameterization of Eq. (p6D, and the relations pi = 0, 
^ 0, and 2P ■ Pc > 0. Applying {d/dkf) to L(/cij,, /c2j,), and letting ki^ [i = 1,2) go 
to zero, the first term in Eq. (^) develops a double pole, while the second term remains a 
single pole. 



_d_ 

dk{ 



L{ki^ = 0, fcgj, = 0) = 7 ■ Pc 



d 



9s 



2P-Pc, 



1 



Xi — X2 + it) 

1 



[xi - X2 + le) 



^^pL{ki^ = 0, k2j, = 0) . 



(52) 



Since we keep only contributions proportional to {d/dx)Tp^\x, x), we neglect the single-pole 
term in Eq. (|5^) in the following discussion, and use 



d 



L{ki^ = 0, k2j, = 0) 



d 



dk{ 

1-Pc 



L{ki^ = 0, k2j, = 0) 



2P-pc 



[Xi - X2+ it) 



(53) 



Similarly, for the diagram at the right in Fig. |TT]b, we have 

1 ■ {Pc + h-ki)'^ ■ P 



R{ki^,k2j,) = Qs 



(Pc + k2- ki) 



It 



1-Pc 



X2- xi+ xo{k2T,kij,) - it 
■ P-f ■ {k2j, - ki^) 



1-Pc 



2P ■ Pc [x2 - xi+ Xo{k2j., ki^) - it] 
where Xq is defined in Eq. (^TJ). Taking the derivative with respect to fcf, we have 



(54) 



d 

^i?(/Ciy = 0, /C2r = 0) 



d 

■^i?(fcl^ = 0,^21- = 0) 



7 - Pc 



2P-Pc 



X2 — Xi — It) 



(55) 



Eqs. (|53D and ( ^51) show that the double-pole contributions from diagrams in Fig. |]b satisfy 
Eq. (^0|). Keeping only these double-pole terms, as in Eq. (42a), we now have 



OK 



(Xi,X2,x',]9c) ~ ^ H2^2{^l,X2,x',Pc) 



2tx 



i^Pc,) 



dH 
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Dr I / \ 9s ttR I I 

-p-\X\^ X2-, X ,Pc) ^ — -"2^21^1; ^2, X ,Pc 



Xi - X2+ it)'^ 

1 



1 



''' ' X2 — X\— itY 



2P-pc 
( 1 



2P-Pc 



(56a) 
, (56b) 
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where Dl and denote the double-pole contributions from the left and right diagram 
in Fig. respectively. In Eq. (|56|) , if|l^2(^i5 ^2, a^',Pc) and if^2(^i) ^2, 2;', Pc) are 2^2 
partonic parts corresponding to the left and right diagrams shown in Fig. |12[ They have 
the limits 



1 



-^2^2(^1) ^2; 3:', Pc) 2:1^x2 — H2^2{X2,x' ,P. 



(57a) 
(57b) 



X2 

Xi 

where H2^2{xi,x',pc) with i = 1,2 are the same as in Eq. (^2]). 

Recalling that Tp is real, it is evident from Eq. (29) that we need the imaginary part 
of {d/dk2)Hj:,^ and {d/dk2)H£,^ in order to get a real contribution to the spin- dependent 
cross section. For double pole terms like those in Eq. ( p6a| ) and ( ^6b|) , the imaginary part 
is given by 



dxi 



1 



xi - X2 + ie) 



F{xi,X2) = / dxi [— ivr 5(xi — X2)] 



d 

dxi 



F{xi,X2) 



(58) 



for any smooth function F{xi,X2)- Using Eq. (0), we have following relation, 



dx,dx2T^/\x^,X2 



d 

-Qj^iHDA^U^2,x',Pc)+HD^{xi,X2,x',Pc)) 
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dx2 
+ J dxi 



d 



dxi 
d 



H^^^{Xi,X2, x',pc) T^^\xi,X2) 



8X2 



H^__^^{Xi,X2, x',pc) T^^\xi, X2 



^srinn fdx a . (v^) 



d 



J Xl=X2 



-I X2=X\ . 



(59) 



In deriving Eq. (p^), we have used the symmetry property T^\xi,X2) = T^^'*(x2, Xi), Eq. 

(H), El 



d_ 

dx 



T^/^\x,x] 



d 
dxi 



{tZ\xi,x) 



J x\=x 



_d_ 

dxi 



(T^/J{x,X2) 



(60a) 



J X2=X 



and Eq. (0). In addition, we have used the approximation 
d 



dxi 



[H^_^^{Xi,X2, x',Pc) TP{Xi,X2] 



J Xl=X2 



— H2^2{x2,x',Pc) 
X2 



d 

dxi 



(61) 



Xl=X2 



demanding as usual a derivative of Tp 

Substituting Eq. ( p9D into the cross section Eq. (^9]), we obtain the leading double-pole 
contributions from the diagrams shown in Fig. 



21 



1 dz 



a c ^mm 



1 c/x 



-T 



G{x') 



^min XjS* ~\~ U I Z 

.— (tZ\x,x) 



dx' ( , 

— - X 



X 



—xT I z 
xS + U/z^ 



(62) 



where T = —2P-i is defined in Eq. (H), and the partonic hard part i, u) is normahzed 



to have 



H. 



D 

ag~*c 



(63) 



with H^g_^^{s, i, it) the same partonic hard part derived from the contribution of the deriva- 



tive of the (5-f unction, and given in Eq. (^). 

In addition to the diagrams in Fig. ^d, there is another type of diagram with final-state 
interactions, as shown in Fig. In this case the final-state interactions taken place on 
an unobserved final-state parton. The diagrams on the left and right are the same, except 
for the final-state propagator and the argument of the phase space 5-function. The total 
partonic contribution from these two diagrams can be expressed as 



Hc^ (fci, k2, x', pc) + Hcj^ih, k2, x', p, 

1 .,,o. 1 



Lf + ie 



m) + 



Ll - It 



F{ki,k2,x',pc) 



(64) 



where the momenta Li and L2 are defined in Eq. (^). The function F{ki,k2,x',Pc) repre- 
sents the common factor of two diagrams in Fig. it has the symmetry property 



F{ki, k2, x', Pc) = F{k2, ki, x', p^) . 



(65) 



/^From Eq. (|5^, combining the symmetry properties of Eqs. ( P^D and (p5D, we readily show 
that the leading contribution of the diagrams in Fig. ^c to the spin- dependent cross section 
(or Eq. (PUp) vanishes. 



C. Quark-Quark and Quark- Ant iquark Subprocesses 



In this subsection, we present the leading contributions to the spin-dependent cross 
section from quark-quark and quark-antiquark subprocesses. 

Based on the same arguments following Eq. (|37|), the leading contributions from diagrams 
with initial-state interactions, shown in Fig. |T3|a, come only from the derivative of the phase 
space 5-function. By analogy to Eq. (^T]) , we obtain 
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STinn 
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—X 



9 / (^>a;) 



dx 



X 




where the partonic hard part, H[g^^{s,t,u) is given by the 2 — > 2 quark-quark (quark- 
antiquark) diagrams shown in Fig. Compared to the spin-averaged case, H^g^^(s,i,u) 
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plays the same role as aaq^c in Eq- 
diagrams needed to calculate daq- 



(p^). In fact, H^g_^^{s, t, u) is given by the same Feynman 
^c, but, with different color factors, C^, due to the extra 
initial-state interactions. Similarly to C^, is given by the color structures of the diagrams 
shown in Fig. p!3|a, contracted with [2/(A^2 — where B and ij are color indices for 

the gluon and quarks from the polarized hadron, respectively. 

Contributions from the derivatives of the phase space 5-functions of the diagrams with 



final-state interactions shown in Fig. 13b are given by 



dH 



ai ^ dz ^ ^ dx 1 

^ a C ■^niin ^ '^^min ^ ~\~ U j Z 



sj'inn 



x'S + T/z^ 



Y,q{x') 




—xT/z 
xS + U/z^ 



-4 



(67) 



where superscript F represents the final-state interactions. The partonic hard-scattering 
function, H^g^^{s,t,u), has the same functional form as Hl^^_^^{s,t,u) in Eq. (|6^), with a 
different color factor C^, because of the final-state interactions. As with C^, is given by 
the color structure of the diagrams shown in Fig. p!^, contracted with [2/(A^2 — 

In addition to the contributions given in Eq. ( ]H7] ) from the derivative of the 5-function, 



the diagrams in Fig. |l^b also have leading contributions from double-pole terms. Just as 
for the contributions from the quark-gluon subprocesses, given in Eq. (|62|), the quark-quark 
and quark-antiquark double-pole contributions take the form 



dH 
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— Dc^^[z 



a,C ^mm 



Z^ 



STinn 



in X xS + U/z 



dx' 
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—xT/z 
xS + U/z^ 



(68) 



where the hard-scattering function found from the double pole, H^^ ^^(s,t,u) is equal to 



Haq- 



(s, i, u) in Eq. 



D. Calculation of the Partonic Hard Scattering Functions 



In Eqs. (^Tf), (^), (|6^), (|66D, ( |67D and (^), we have presented factorized expressions for 
leading contributions to the spin-dependent cross section, EidAa{sT)/d^i, for quark-gluon, 
quark-quark and quark-antiquark subprocesses. To complete our derivation of the spin- 
dependent cross section, in this subsection we outline the calculation of the partonic hard 
scattering functions H^^^^, H^^^^, H^^^^ and H^^^^. We recall that the subscripts I and 
F refer to initial- and final-state interactions, respectively. The other two hard scattering 
functions, associated with derivatives on final-state propagators only, Hag-*c ^"^^ ^aq^c^ 
equal to Hag^c ^aq^c^ respectively. 

For the quark-gluon subprocesses, the partonic hard scattering functions H^^^^ and 



Han^c given by the same quark-gluon 2 — >■ 2 Feynman diagrams as shown in Fig. [15 



which are actually the same diagrams contributing to the spin-averaged partonic part, aag^c, 
in Eq. ( |35a|) . Incoming quark lines are contracted by (1/2)7- {xP)^ and incoming gluon lines 
are contracted by (1/2) (—(^q,^). 
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Let Cg, Cg, and Cg be the color factors for processes that are spin-averaged, spin- 
dependent with an initial-state interaction, and spin-dependent with a final-state interaction, 
respectively. The factor Cg for each diagram shown in Fig. |15| is simply the standard color 
factor for that diagram, with an average over initial-state quark and gluon color. 

Each Cg is given by the color factor of the diagram with one extra initial-state three- 
gluon vertex. An example is shown in Fig. |l6|a. The color of the incoming gluon from the 
unpolarized hadron is averaged, and the colors of the incoming quarks and the extra gluon 
from the polarized hadron are contracted with [—2i/{N'^ — l)]{t^)ij, as explained in the text 
following Eq. 

Finally, the Cg are the color factors of the same 2^2 diagrams with an extra final- 
state quark-gluon interaction, illustrated by the diagram shown in Fig. |l6lb. Similarly to 
Cg, the color of the incoming gluon from the unpolarized hadron is averaged, and the colors 
of the incoming quarks and the extra gluon from the polarized hadron are contracted with 
[2/{N'^ — l)](t^)ij, as mentioned in the text after Eq. (^Uj). Our results for all these color 
factors are collected in Table |. 

For quark-quark (or quark-antiquark) subprocesses, the partonic hard scattering func- 
tions, -fffg^c -^og^c given by the same quark-quark (or quark-antiquark) 2 — 2 
Feynman diagrams as shown in Fig. |T^, which are the same diagrams contributing to the 
spin-averaged partonic cross section, aaq^c in Eq- (|35b| ). Incoming quark lines from the 
polarized hadron are contracted by (1/2)7 ' {xP)^ aiid incoming quark (or antiquark) lines 
from the unpolarized hadron are contracted by (1/2)7 ' {x'P'). 

As with the quark-gluon subprocesses, Cq, C^, and C^ are respectively the color factors 
for subprocesses that are spin-averaged, spin- dependent with an initial-state interaction, and 
spin-dependent with a final-state interaction. The Cq for the individual diagrams shown in 
Fig. are the color factors for each diagram, with a standard average over initial-state 
quark (or antiquark) color. The C^'s are found by including an extra initial-state three- 
gluon interaction in the 2^2 process, (for example. Fig. [!%.) averaging the color of the 
quark (or antiquark) from the unpolarized hadron, and contracting the colors of the incoming 
quarks and the extra gluon from the polarized hadron with [2/(A^^ — l)](t^)jj, as mentioned 
following Eq. (|66D. The C^ are found from the same 2^2 diagrams, now with one extra 



final-state quark-gluon interaction (illustrated by the diagram shown in Fig. |T8|b). In exactly 
the same fashion as for C^, the colors from the unpolarized hadron are averaged, and the 
colors from the polarized hadron are contracted with [2/(A^^ — l)](t^)jj (as mentioned in 
connection with Eq. (|B7D). 

Our results for the quark-quark and quark-antiquark color factors are summarized in 
Table |T[ Notice the sign difference for the coefficient of AN in the color factor (A^^ ± 
AN — 4)/ (32A^), between graphs related by reversing the arrow of a quark or antiquark line. 
These will give slight differences to the asymmetries in proton (|) -proton compared with 
antiproton(|)-proton collisions. 

/^From Table | and Table |I|, we can construct all the necessary partonic hard scattering 
functions. For the spin-averaged cross section, the hard-scattering function for the quark- 
gluon subprocess, dag^c in Eq. ( p5a|) , is found by combining the entries in the columns of 
Partonic Parts and Cg in Table |I[ For the quark-quark (or antiquark) subprocesses, CFaq^c 
in Eq. (|35b| ) is found by combining entries from the columns of Partonic Parts with Cq 
in Table fl[ For the spin-dependent cross section, the twist-three partonic hard scattering 
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function H^g^^ is found by combining entries in the columns of Partonic Parts and in 
Table |. In the same way, one can read off other partonic hard scattering functions, H^g^^, 
Hag^c and H^^^^ from Table | and Table |I|. 

IV. NUMERICAL RESULTS FOR SINGLE TRANSVERSE-SPIN 

ASYMMETRIES 

Having derived expressions for the single transverse-spin asymmetries in previous section, 
we are now ready to develop numerical estimates of Aj\f for inclusive single pion production. 

A. Model for the Twist-3 Distribution: T^\x,x) 

The application of perturbative QCD to observables involving hadrons in the initial 
state relies on factorization theorems ||31[ and on the universality of the nonperturbative, 
long-distance distributions. For the single transverse-spin asymmetries discussed in this 
paper, a test of the perturbative formalism requires in principle an independent extraction 
of the spin-dependent twist-three distributions, (j)^^/j^{xi,X2) introduced in Eq. (0). As we 
have observed, there are a variety of twist-three distributions, dependent in general on a 
pair of momentum fractions. It would requires extensive measurements to pin down all of 
these functions. However, for single-spin asymmetries in the forward region, we have argued 



above, and in Ref. |fOl, that the dominant contribution may depend primarily on only a 



single twist-three distribution, T^\x,x), at equal values of its two arguments. Assuming 

this to be the case, it could be possible to infer the form of T^\x, x) from single transverse- 
spin asymmetries in tt"*" and/or vr" production, and then use it to predict asymmetries in 
the production of 7r°, direct photon or other particles, at least approximately. 

In order to compare our calculated asymmetries to the existing data, we need to assume 
an initial functional form for the twist-3 distribution, T^'' (x, x). To help motivate our model, 
we compare the operator definition of Tjy\x,x) with that of a twist-2 quark distribution 
qa{x) of flavor a. From Eq. (^), we have 



J dy^ e^--"^ F^+iy^)^ My-)\P, ^t) , (69) 

where subscript a is quark flavor. Correspondingly, from Ref. |39] we have for the quark 
distribution 

qaix) = J ^e^^^"^- {p\Moh+My-)\P) ■ (70) 



As above, we suppress ordered exponentials of the gauge field. Comparing Eqs. ( |69D with 
([fO|), the operator defining the same as for the spin-averaged quark distribu- 

tion, except for the term in the square brackets. This factor, however, does not introduce 
explicit x-dependence (or ^/-dependence in coordinate space). Based on this similarity of the 



25 



operators, we model the twist-3 distribution with the following functional form, inspired by 
the quark distributions themselves. 



Tj^^ {X, X) = KaX Qa (x) , (71) 

where A (with dimensions of energy) is a normalization constant, which will be fixed by the 
data; and where Ka = ±1,0, depending on flavor a. Note that we propose the relation Eq. 
(fflp only for relatively large x, where the correlations of quarks with the gluon field may be 
simplified. This restriction limits somewhat the utility of low moments of Tp in estimates 



of its magnitude [|T^] 



For the parameters Ha in Eq- (|7l|) , we shall see that the data suggest the choices 



Ku = +1 and — = — 1 (proton) , 

= —I and — = — 1 (antiproton) , (72) 

where the second line follows from the first by using charge conjugation invariance in Tp. In 
the valence quark approximation, discussed in the previous sections, we further assume that 
Ks = 0. Of course, Eq. (|7TD is simply a model, and the true functional form of the twist-three 



distribution T^\x,x) should be determined by detailed comparison with experiment. The 
purpose of our model is to have a functional form that we can use to begin such a comparison 
with the important, but still limited, data that are available. 



B. Single Transverse-spin Asymmetries in Pion Production 

Single transverse-spin asymmetries for pions were been measured at Fermilab by the 
E704 Collaboration with 200 GeV polarized proton and antiproton beams on an unpolarized 
proton target 0]. In this subsection, we use the Fermilab data to estimate the value of A, 
in Eq. (j7I|), and check the consistency of our model. 



1. Absolute Sign of the Single Transverse- Spin Asymmetry 



In order to compare the experimental data on the asymmetries. An, with our calculations 



in Sec. |T|, we need to fix the absolute sign of A at. 

According to Ref. 0], positive values of An correspond to larger cross sections for pro- 
duction of 7r° to the beam's left when the beam particle spin is vertically upward, as sketched 



in Fig. 0. We choose our coordinate system such that the beam direction is along the z-axis, 
and the direction of the beam spin is along the shown in Fig. [T^. Consequently, 

the experimental beam's left corresponds to the — ^/-direction in our coordinate system, and 



(A;v)™ > 



/cxp 



> . 



(73) 



Eq. (|73| ) fixes the absolute sign oi An presented in Sec. |T|, and dictates our choice k^j = +1 
in Eq. (|7l). 
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2. Leading Single Transverse- Spin Asymmetry [{d/dx)T^\x,x) only] 



As explained in Sec. |T|, we are interested in An in the forward region, where it is 
largest experimentally. In deriving Eqs. (^7D, (^), (p^, ( pTD and (p^), we kept only 

contributions from the terms discussed in items (1) and (2) of Sec. [II D| , because those 
discussed in items (3) and (4) lack a derivative on the twist-three distribution. To be 
consistent with our approximation, we rewrite the contributions in Eqs. (^Tj), (^91), (|62D, 
(|66|), (p^) and (|68D in terms of an explicit factor of {d/dx)Tjy\x,x), neglecting derivatives 
of other factors. Combining all leading contributions to the spin-dependent cross section, in 
a manner similar to the spin-averaged cross section in Eq. (B3), we obtain 
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Gix')Aaag^, +Y,q{x')Aa. 



aq^c 



where J2a runs over up and down valence quarks. The integration limits in Eq. (^4|) are the 
same as those defined in Eq. (pSj). The spin-dependent partonic cross sections, A&ag^c and 
Aaag^c are given by 



ag^c 



Hi S ,A , ZJF 



ag^C 



Aa, 



aq^c 



Hi(sXu) + Hi 



u 



t 



u 



H. 



D 



S, t, u) 
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(75a) 
(75b) 



where the minus sign is from e^^^nri _ _ ^esT-nn^ where all the partonic hard-scattering 
functions have been given in Sec. fT|. In deriving Eq. (|75D , {x'S + T/z)/{—T/z) = u/t was 
used. From the information given in Table | and Table we find the following explicit 
expression for Aaag-,c taking = 3, 
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(76a) 



For quark-quark (or antiquark) scattering, the color factors of individual subprocess depend 
on quark or antiquark, as shown in Table ||. For parton a a quark (corresponding to a 
polarized proton beam), we have 



A&qq'^q 



4 ( + 



9 



t2 



64 8 V t , 
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For a polarized antiproton beam, 
ble 0. 



similar formulas for Aa^fc^c can be derived from Ta- 
we see that the underlying partonic cross sections for 



From Eqs. (|33D and 

spin-dependent and spin-averaged cases are similar, other than the factors in square brack- 
ets. 



3. Comparison with the Fermilab Data 



Because of limited phase space, most of the Fermilab data in Ref. 0] were collected at 
relatively small values of transverse momenta, ranging up to 4 GeV for tt^ in the central 
region (where Aj^ is small), and up to only 1.5 GeV for tt"^, tt^ in the forward region, where 
Ajv is large. In general, a transverse momentum of even 2 GeV is considered too small 
to apply perturbative QCD reliably to single-particle inclusive cross sections, because of 
their steep dependence on Et- This strong dependence makes the cross sections sensitive to 
higher-twist effects not associated directly with spin, such as intrinsic transverse momentum, 
hadronic scales, and, of course, yet higher powers in l/ix- One consequence of these effects is 
to regularize the cross section at ix = 0. For the asymmetry, however, the strongest power 
dependence on 1/£t cancels in the ratio of the spin-dependent and spin-averaged cross 
sections, leaving at most A/£t in Atv- In fact, as we will show below. An does not behave 
numerically even as steeply as 1/£t in most of the range where the data were collected. 
This suggests that our calculation for is perturbatively stable and may be meaningfully 
compared with the data. 

For simplicity in our numerical estimates, we employed the following simple parametriza- 



tions, without scaling violation, for twist-two parton distributions ||T0 
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Here, B{x,y) is the beta function. For pion fragmentation functions, we rely on Ref. 



Using the simphfied parton distributions of Eqs. (|77D in the spin-averaged cross section, and 
in the model for the twist-three distribution given by Eq. ([7ID, we evaluated as the ratio 
of the spin-dependent cross section in Eq. (|7^ to the spin-averaged cross section, Eq. (p3|). 

In Fig. |2y, along with experimental data from Ref. Q], we have plotted our calculated An 
for 77+ and vr" production in the scattering of a polarized antiproton beam on an unpolarized 
proton target. Similarly, in Fig. we plot the asymmetries with a polarized proton beam. 
In Fig. we compare theory and experiment in the asymmetries for vr^ production with 
a polarized antiproton beam and a polarized proton beam. The data presented in Figs. 
pT| and ^ are averaged over the range of transverse momenta, up to 1.5 GeV. All of the 
calculations in these figures, however, were carried out at ~ 4 GeV, with a normalization 
constant A = 0.080 GeV, adjusted to give a rough match to the data. [| We will come back 
to the choice of £t in a moment. This limitation notwithstanding, fixing the single overall 
normalization constant. A, is enough to give theoretical predictions that are consistent with 
the shapes and relative signs and normalizations of all the experimental data. 

Now let us consider to the question of how to best to choose ix for the comparison with 
the data. Given the naive expectation that An ~ the extracted value of A might 

be expected to depend strongly on the value of It at which we evaluate the asymmetries. 
Surprisingly, however, the perturbative prediction for the asymmetries in this momentum 



region is not very sensitive the precise value of Et- Thus, in Figs. 23, E4l and 25, we present 



the same asymmetries as in the foregoing three figures, now evaluated at ix = 1-5 GeV. 
For this value, we find a good match to the data by choosing A = 0.070 GeV, not too 
different from the value found at £t = 4 GeV. Clearly, the normalizations and shapes of 
the asymmetries at £r = 4 GeV and = 1-5 GeV are very similar, with an only slightly 
different normalization factor. We consider this stability very encouraging. Such consistency 
is strong evidence that the twist-three formalism of perturbative QCD can be applied to 
single transverse-spin asymmetries at moderate transverse momenta. We will give a further 
discussion of this point in the next section. 

We close this section with a few comments on the consequences of our model of Tp 
(Eqs. ([7I| ) and (|72|), with A ~ 0.080 GeV) for single-spin asymmetries in direct photon 
production. Compared to the ansatz for Tp proposed in Ref. fl^, the two main differences 
are, first, the relative minus sign between the down and up quark matrix elements, and, 
second, a decrease in the overall normalization A, below 100 MeV. Both of these features 
are suggested by comparison to the data for pion production, which is only now possible. 
The effects of the both changes would be to reduce the cross section estimates given in Ref. 
n| , although the second is more important than the first, because the down quark's charge 



is small. In any case, the data which limit the direct photon asymmetry is at low xp, 
where either model predicts a small effect. 



For the purpose of this comparison, we neglect correlations between xp and ip in the data. 



29 



V. SUMMARY AND DISCUSSION 



In this section, we summarize and interpret the main features of our results and provide 
a few thoughts on future development on this subject. 

We have presented a calculation of single transverse-spin asymmetries, A^r, for hadronic 
pion production at large xp- This calculation was based on a "valence quark-soft gluon" 
approximation. In this approximation, we kept only those contributions to An proportional 
to the derivative of the twist-3 quark gluon correlation function, {d / dx)T^\x , x) , where 
a denotes a valence quark flavor, and where the equal arguments in Tp imply zero gluon 
momentum fraction. Our results for spin-dependent single-spin cross sections are given in 
Eq. ([7^). The ratio of the spin-dependent cross section in Eq. ([7^ ) and the spin-averaged 
cross section in Eq. (^) defines An for hadronic pion production. The spin- dependent 
cross section, Eq. ( |7^ has two types of contributions: quark-gluon and quark-quark (or 
antiquark), which are given by Adag^c in Eq. ( [76 a|) and A&ab^c in Eq- (|76b| ) , respectively. 



All of these calculations are strictly leading order; we anticipate that a large part of higher 
order corrections will cancel in the asymmetry. Our model for the twist-3 matrix element 
Tp is given in Eqs. (^) and ([72[). We have not investigated the evolution properties of these 
matrix elements here. We expect this to be an interesting subject, but we do not anticipate 
that evolution will require qualitative changes in our conclusions. 

Single transverse-spin asymmetries are a twist-three effect in QCD perturbation the- 
ory. After taking the ratio of Eqs. (|7^ and (|33D , the asymmetry has following schematic 
dependence on kinematic variables in the large xp region. 



A 



N 



It 



-u) 



\ — Xp 



(78) 



where the invariants, U and T are defined in Eq. (||). In Eq. ([78|) , the prefactor ix/i—U) 
comes directly from the factor e^''^""/(— -u) in the spin-dependent cross section in Eq. ([741). 
The combination [1 -|- 0{U/T)] is left over from the partonic cross sections in Eq. ([f6[), 
after the cancelation of the dominant dependence in the ratio. The normalization 
parameter A comes from our model of the twist-three correlation functions, Tj^\x,x) in 
Eq. ([TTp. Finally, the factor 1/(1 — xp) for xp large is associated with {d/dx)T^\x,x) 
in Eq. ([7^). The approximate 1/(1 — xp) behavior in the ratio of the derivative of the 
twist-three correlation function to the corresponding twist-two parton distribution is the 
dominant feature of the twist-three asymmetry, and is responsible for the observed growth 
of An in the large-x^? region. 

The factors ix/ {~U) and ip/i—T) in Eq. ([f8[ ) reflect the twist-three nature of the asym- 
metry. An- Combining Eqs. (^ and (H), we express the invariants, U and T, in terms ofxp 
and ip, 



U = - 
T = - 



Xp ~\~ X'p ~\~ Xp 



™2 1^ ™2 



Xp 



(79a) 
(79b) 



When Xp = 0, both U and T are equal to ipVS. /,From Eq. (^), we conclude that the 
asymmetry a.t xp = should have a very mild, probably linear dependence on the pion's 
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transverse momentum (A^ must vanish at ix = 0). Our analytical results in Eq. ([f4D, 
however, are not accurate for the asymmetry near xp = 0, because of the large xp approx- 
imation used in our calculations. But, from the general structure of the asymmetry, we 
believe that weak transverse momentum dependence at = for should be a more 
general conclusion. 

U xp ^ xt, the invariants U and T in Eq. (^) have the following approximate depen- 
dence on ip and Xp, 

U -Xp S , (80a) 
T ^ . (80b) 

Xp 

Consequently, in the large xp region, the asymmetry, Aj^, will have two typical contribu- 
tions, X/ip and Xip/S, respectively. If the X/ip contribution dominates, perturbative QCD 
calculations of the asymmetry may be relatively sensitive to nonperturbative effects, because 
of its singular behavior at ip = 0. On the other hand, QCD perturbation theory may pro- 
vide a reliable calculation of the asymmetries when the Xip/S term is relatively important. 



In Fig. ^ we plot the transverse momentum dependence of the asymmetry at xp = 0.4, 
where most data were collected. The asymmetries for both tt~^ and vr" have a quite weak 
dependence on pion's transverse momentum for ip > 2 GeV. This suggests that perturbative 
calculations for the asymmetries are reliable for a wide range of the experimental kinematics. 

The remarkable feature of mild transverse momentum dependence, shown in Fig. |2^, 
can be easily traced to Eq. (^). For the quark-gluon subprocess, once the dominant 
dependence has canceled in the asymmetry, the coefficient of u/t is much smaller than 
the corresponding constant term. Similarly, for quark-quark and antiquark subprocess, the 
coefficient of u/t is also much smaller than the constant term, except for terms proportional 
to l/it^ and The latter, however, are suppressed by t^/s^ relative to the leading terms 
in the forward region. In summary, the small coefficients for u/t terms assure that the 
ip/i—T) dependence in Eq. ( [THD does not dominate the ip/{—U) dependence. We verify 
this conclusion by plotting the fractional contributions to the vr^ asymmetry from l/{—U) 
term and l/(— T) term, respectively, as a function of pion's transverse momentum in Fig. |2^. 
It is evident that contribution from l/{—U) term is comparable with l/(— T) term for the 
region of our interest. 

If Xp — > 1, or U/T x"]^ S/ip ^ 1, the asymmetry will be eventually dominated by the 
X/ip terms. Therefore, the asymmetry will scale with 1/ Ipin this region. For the kinematics 
of the Fermilab data, this scaling region is not yet reached. In Fig. the steep increase 
of the asymmetries for ip < 2 GeV indicates the dominance of the X/ip contribution, and 
probably signals that the perturbative calculations are relatively less reliable if ip is much 
less than 2 GeV. The slight increase when ip ^ 6 GeV signals an effect of the edge of phase 
space. Clearly, the high energies of the polarized RHIC proton beam would make it possible 
to check these predictions. In Fig. ^ we show the ^^-dependence of A]\r for xp = 0.4 at 
y/S = 200 GeV. Compared to Fig. |2| at Fermilab energies, the ip/U term is relatively 
suppressed, and the model predicts a steeper ^^-dependence and, in general, a smaller, but 
still substantial, asymmetry. Fig. |2D| shows the asymmetry as a function of xp at ip = 4 
GeV. These are examples only; the model can be used to predict A^- over any kinematic 
range that is experimentally convenient, so long as it is in the forward region. 
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In summary, we have calculated the single transverse-spin asymmetry for hadronic pion 
production in perturbative QCD. With only one normalization parameter A and a relative 
sign of polarized twist-3 valence quark distributions, our numerical results are consistent 
with Fermilab data on the asymmetry for both the sign and shape, as well as relative nor- 
malizations. In addition, we demonstrated that perturbative calculation of the asymmetries 
is applicable even for pion momenta as small as a few GeV. This conclusion is very encour- 
aging for future applications of perturbative QCD beyond the leading twist. Our method 
can be easily generalized to calculate the single transverse-spin asymmetries for inclusive 
production of other particles. The planned polarized beam at RHIC affords an exciting 
opportunity to test these, and related ideas on the spin structure of the nucleon. 
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APPENDIX 

In this appendix, we identify twist-3 distributions and fragmentation functions that can 
contribute to the sums in the collinear expansion, Eq. ([T7|). The factorization in ( p!?] ) enables 
us to apply parity and time-reversal (PT) invariance to hadron-hadron scattering in a manner 
similar to their classic application to inclusive DIS, reviewed in Sec. |11 A|. Thus, it will be 



natural to study the symmetry properties of possible matrix elements. 

We will identify terms of the type discussed in connection with Eq. (^), that is, with 
integrals over two quark momentum fractions, Xi and X2- Furthermore, we require that Xi 
be set equal to X2 by a "gluonic" pole at Xi = X2 in the hard scattering (see Eq. (^7D), in 
accordance with our valence quark-soft gluon approximation. Let us concentrate first on 
parton distributions, and return at the end to fragmentation functions. 



Twist-3 Distributions 



As mentioned in Sec. |11 C| , the derivation of terms in Eq. (0) involving quarks requires a 



Fierz projection of the Dirac indices linking the distribution or fragmentation function and 
the hard scattering. A schematic illustration was given in Fig. |^. The collinear expansion 
then isolates twist-three fermion matrix elements with two quark fields and either a covariant 
derivative or a field strength. It will be convenient to start by discussing the expectation 
values of combinations of these fields in position space. We thus introduce 

D^Avi, 1/2, s) = (P, s\ m r D\y2) i^iyi) \P, s) (81) 
F^iVu 2/2, s) = {P, s\ m r n.r^^iy^) ^yi) \P, s) , (82) 

with r a Dirac matrix. We define = id^ —gA"^, and we adopt the kinematics and notation 
of Sec. |1I A| ; in particular, is defined in Eq. (^. In these matrix elements, the index i 
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is assumed to be transverse. This alone is enough to make the matrix element twist-3; the 
Dirac projection must not raise the twist further. The relevant terms in the Fierz projection 



between the distribution for a hadron of momentum P'^ = and the hard scattering 

are then given by 

^aa'^bb' = ^{l ■ n)ab{l ■ n)b'a' + " ^75)afe(757 " ^)b'a' 

+^E(M^UMO,,„, + ---, (83) 

where omitted terms raise the twist, and where we define 

(na)'^ = n^a^"^ . (84) 



For an opposite-moving hadron, with momentum P'^ = we exchange the roles of 

and n^. The matrices above have the properties 

r = 7° rt 7° , (85) 

T = Sr (Tr*T)^ , (86) 

with 5r = ±1, where T = i'-f^'-f^ = is a time-reversal matrix that acts as 

T(7^)*r = 7^. (87) 

Specifically, for the vector, axial-vector and tensor cases we have 

n ■ 7 : Sn-^ = 1 (88) 
n ■ 775 : 5„.775 = -1 (89) 
{nay-. 5(^nar = -l. (90) 

^From the expectation values and Dp we define parton distributions by Fourier transforms 
with respect to light-cone momenta, and if desired transverse momenta as well, 

t?^\x,,k„X2M,s) = J dy,dy2e"'-y^+'^'^-'^>ywUyi,y2,s), (91) 
tPix,M,X2M,s) = J dy,dy2 e''''-y^+'^>''-'^>y'Fliy,,y2,s) , (92) 

where we define dyi = dy~d^y, with y a two-dimensional transverse vector, and ki ■ y = 
XiPy~ ~ k ■ y- In the following, we study constraints on spin-dependence that follow from 
the reality and symmetry properties of these matrix elements in QCD. This will enable us 
to identify the relevant contributions to the sums in Eq. (\V7\j . 

Reality and Symmetry 

The reality properties of the matrix elements (0) and (|82D are conveniently expressed 



as 
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'Vi, 



-^2, S) 



Dr{yi,yi - 1/2, s) 
Fr{yi,yi-y2,s) 



which relate, of course, expectation values with the same spins, 
reversal and parity, on the other hand imply that 

-D^(?/l, y2, s) = 6r D].{yi, yi - y2, -s) , 
Friyu Z/2, s) = -5ri^r(?/i> l/i - l/2> s) , 

in which spins are reversed. 



(93) 
(94) 

Invar iance under time 

(95) 
(96) 



Relations for parton distributions t\P^^ and tp^"** are easy to derive by inserting the reality 



and symmetry relations into the Fourier transforms of Eqs. (pT|) and (|92|) , and changing 
integration variables. Because in this paper we are concentrating on the coUinear expansion, 
with convolutions in light-cone momenta only, we shall suppress transverse momenta in the 
arguments of the distributions, and exhibit only the momentum fraction variables Xi in the 
following formulas. Relations for transverse-momentum distributions are found by simply 
reinserting the kj arguments, alongside the corresponding momentum fractions. With this 
understood, the reality conditions give 



4^^* {xi,X2,s) = t^'"" {X2,XI,S 



{D)i. 



t[f^*(Xi,X2,s) J =t''^^\x2,Xi,s) , 

in which we note that the momentum arguments are exchanged. Referring to Eqs. 
90) above, we see that the even parts of the twist-3 distributions are real, the odd parts 
imaginary. 

Similarly, from PT invariance, we find 



(97) 
(98) 



tr {xi,X2, s) 
i:|f^'(xi,X2, s) 



6r t^j^^\x2,xi, -s) , 
-6r t\^^\x2,xi, -s) 



(99) 
(100) 



with 6r defined in Eq. (^). Note the extra minus sign in the second case, which reflects the 
PT properties of the field strength tensor. 

/^From Eqs. (|97| - |100|) , we can derive the constraints on the spin-averaged. 



{tp^'){Xi,X2) = ^ [t'^'\xi,X2,s) +t''^'\xi,X2,-s) 



and spin-dependent 



AtP(Xi,X2) 



1 

2 



t[?^*(Xi, X2, s) - t''-^'\xi,X2, -S) 



(0)i, 



(101) 



(102) 



distributions for each choice of operator O = D,F and Dirac structure F. Specifically, the 
spin-dependent distributions Atl^^% At^^^*^ and At|^^*^j and the spin-averaged distributions 

(^1^75); {^[nly) (^r^7*) imaginary and vanish at Xi = X2. They therefore cannot be 
associated with gluon poles in Eq. (|T7|) , and are nonleading in the valence quark-soft gluon 
approximation introduced in Sec. |IIB| . 
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Leading Terms at Twist-3 



The remaining distributions are real and nonzero at Xi = X2 in general. For the first 
sum in Eq. (|17D, we need a real, chiral-even, spin-dependent parton distribution. The only 
one is Atl^^^{xi,X2), which is equal, up to a constant, to Tj^\ Eq. p5|). 



X2j 



(103) 



where the tensor structure follows from parity invariance applied to the matrix element. 

For the second sum in Eq. (|T7|) , we need a chiral-odd spin-averaged distribution, to 
give a nonzero trace in the hard-scattering amplitude when paired with the transversity 



distribution 22 



6q{x) 



dy 
27r 



(P,s|V^(0) -{naYsa^ij{y-n) \P,s) . 



(104) 



Here again there is only a single contribution, (t(^^)j)(a;i, 0:2). Parity invariance implies that 
(^(^])j )(^i'^2) is of the form 



(t[£;,)(xi,X2)=4vr:|Tr(a:i,X2), 
where the scalar distribution T'f^ is defined by analogy to Tp'\ Eq. (^ 



(105) 



Tj^\xi,X2) 



dyidy2 ^.xiP'-y++i{x2-xi)P'-y+ 

4:71 

xi 5:(P', .'I ^(0) (na), n,F^'{yt) ^{yt) \P\ s') . 



(106) 



In these expressions we have taken P'^ in the minus-z direction, in accordance with the 
kinematics of the unpolarized hadron in Sec. |II A| . In Ref. |]10| the possibility of such a term 
was noted. 



Fragmentation at Twist-3 



Turning to the third term in Eq. (|T3), we must deal with twist-3 chiral-even fragmentation 
functions, which are transforms of matrix elements of the general form 



d»(l/i,l/2,^) = E Tr f (0| V^(O) |£,X)(£,X| D\y2) ^(j/i) |0) 



X 



/('^^ (2/1, 1/2,^) = E Tr [ (n,a), (0| ^(0) n,^F^\y2) ^(2/1) |0) 

X 



(107) 



(108) 



with the sum over inclusive final (out) states |X, £), where ^ is the momentum of the observed 
particle. The vector is defined by analogy to in Eq. as a lightlike velocity vector 
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in the direction opposite to ne = £^/io- The trace is over Dirac indices. There is no analog 
of the spin variable in this case, although extensions to production of polarized particles 
should be straightforward. We have used the constraints of parity in forming scalar 
fragmentation functions, depending on two momentum fractions. In momentum space they 
are 

S'^\z,,Z2) = [ '^yi^^-^^-ny./z.-U-ny.H/z.-l/z,) d(^)^y^^y^J) (109) 

J An 

f^-\z^,Z2) = I ^yi^e-'^-nyi/^i-^^-ny,il/.,^l/z,) fi<^)(^y^^y^J) . (HQ) 

J Atc 

The constraints of reality are different for these fragmentation functions than for the 
distributions, because the sums over states in Eqs ( |107D and ( |108D are incomplete. We find 
that 

d^''^*i-yi,-y2j) = E Tr [ (nea), (0| ^(0) D\y, - y,) %X){i,X\ ij{y,) |0) ] 

X 

(111) 

/('^)*(-yi, -y2,i) = E Tr [ {nea)i (0| ^(0) n,F^\yi - y^) |^,X)(£,X| i;{y,) |0) ] . 

(112) 



As Collins has emphasized |[T7], time-reversal does not constrain fragmentation functions 



in the same manner as distributions, because T reverses the roles of in and out states. To 
the extent that a sum over in and out states is the same in these functions, symmetry 
under PT would imply that d^"^ is purely imaginary, while /*^°"^ is real. These properties 
can, however, be modified by phases associated with final state interactions. Indeed, this 
is the mechanism by which Artru et al. [|r^ derive single-spin asymmetries starting from a 
model for fragmentation functions with intrinsic transverse momenta. Such functions can be 
thought of as extensions of d^"\ finite distances from the light cone. Following the procedure 
of Sec. ^ above, we can derive hard-scattering coefficients for either function. We reserve 
this for future investigation. 
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FIGURES 




FIG. 1. Inclusive lepton-hadron deep- inelastic scattering, with the target hadron polarized 
transversely. 




+ C.C. 



FIG. 2. Quark-quark scattering diagrams that give a nonvanishing single transverse-spin asym- 
metry in large-pT reactions 




FIG. 3. Sketch of single pion production in spin-averaged hadron-hadron collisions. 
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FIG. 4. Factorization of a typical forward scattering amplitude contributing to the 
spin-dependent cross section for hadronic pion production: (a) with chiral-even three-parton matrix 
element, (b) with chiral-odd transversity. 




FIG. 5. General Feynman diagrams contributing to the partonic parts H in Eq. ([T 
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FIG. 6. Three classes of quark-gluon diagrams contributing to the spin-dependent cross section 
A(t(st): (a) diagrams with an initial-state pole, (b) and (c) diagrams with a final-state pole. 
Symbols B and ij are color indices for the gluon and quarks. 
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FIG. 7. General diagram that gives a leading contribution to Aa{sT): (a) before separation 
of spinor trace and Lorentz indices, (b) leading contribution after separation of spinor trace and 
Lorentz indices. 
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FIG. 8. Factorization of a general diagram contributing to Sa{ki, ^2) of Eq. (^): (a) separation 
of target hadron, (b) separation of final-state pion. 
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(b) 

FIG. 9. Sample leading order Feynman diagrams contributing to the cross section of hadronic 
single pion production. 
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FIG. 10. Two-parton forward scattering amplitude contributing to the partonic hard part 



H2^2{x,x',pc) in Eq. (43). 
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FIG. 11. Sketch for the effective diagrams giving the leading poles in Eq. (^8[): (a) pole to the 
left of the cut; (b) pole to the right of the cut. 




FIG. 12. Effective quark-gluon 2 — > 2 diagrams with the thin line of momentum (x2 — xi)P 
representing momentum flow that is a result of the extra final-state interaction. 




(c) 



FIG. 13. Three classes of quark-quark (or antiquark) diagrams contributing to the 
spin-dependent cross section A.a{sT)'- (a) diagrams with an initial-state pole, (b) and (c) diagrams 
with a final-state pole. Symbols B and ij are color indices for the gluon and quarks. 
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FIG. 14. Effective quark-quark (and antiquark) 2^2 diagrams contributing to the partonic 
hard parts, Haq^c- 
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FIG. 15. All 2^2 quark-gluon diagrams contributing to partonic hard parts, -ff, 
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FIG. 16. Sample diagrams with initial-state and final-state interactions, used to calculate the 
color factors, Cg and Cg in Table |. 
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FIG. 18. Sample diagrams with initial-state and final-state interactions, used to calculate the 
color factors, Cg and Cg in Table |l|. 
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FIG. 19. Sketch for the coordinate system: the polarized beam is along the z-axis and the beam 
particle spin along the x-axis. Positive A]\[ corresponds to an excess of events in the — y-direction. 
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FIG. 20. Single transverse-spin asymmetry as a function of Xi? for 7r+ and vr" production with a 
polarized antiproton beam. Here and in the following five figures, data are from Ref. Q at -v/S* = 20 
GeV and It up to 1.5 GeV. Theory curves are evaluated at transverse momentum /j- = 4 GeV and 
A = 0.080 GeV at the same center-of-mass energy. 
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FIG. 21. Single transverse-spin asymmetry as a function of xp for vr^ and vr^ production 
with a polarized proton beam. Data are from Ref. Q. Theory curves are evaluated at transverse 
momentum = 4 GeV and with A = 0.0.80 GeV. 



46 




FIG. 22. Single transverse-spin asymmetry as a function of xp for vr production with a po- 
larized antiproton and proton beams. Data are from Ref. Q. Theory curves are evaluated at 
transverse momentum /j^ = 4 GeV and with A = 0.080 GeV. 




FIG. 23. Single transverse-spin asymmetry as a function xp for vr^ and vr production with 
a polarized antiproton beam. Data are from Ref. |Q. Theory curves are evaluated at It = 1-5 GeV 
and A = 0.070 GeV. 
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FIG. 24. Single transverse-spin asymmetry as a function xp for vr^ and vr production with 
a polarized proton beam. Data are from Ref. Theory curves are evaluated at px = 1-5 GeV 
and A = 0.070 GeV. 
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FIG. 25. Single transverse-spin asymmetry as a function of x_f for vr production with a polar- 
ized antiproton, along with the same asymmetry obtained with a polarized proton beam. Data are 
from Ref. Q. Theory curves are evaluated at It = 1-5 GeV and A = 0.070 GeV. 
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FIG. 26. Single transverse-spin asymmetry for tt"*" and it production with a polarized proton 
beam as a function of pion's transverse momentum It- Theory curves are evaluated at xp = 0.4, 
VS = 20 GeV and A = 0.080 GeV. 
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FIG. 27. Fractional contribution from !/(—[/) and l/(— T) terms to the single transverse-spin 
asymmetry of tt'^ production as a function of pion transverse momentum. Theory curves are 
evaluated at xp = 0.4, = 20 GeV and A = 0.080 GeV. 
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FIG. 28. Single transverse-spin asymmetry for tt"*", tt and 7r° production with a polarized 
proton beam, as a function of pion transverse momentum It- Theory curves are evaluated at 
XF = 0.4, V5 = 200 GeV and A = 0.0.80 GeV. 
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FIG. 29. Single transverse-spin asymmetry for vr^, vr and vr'' production with a polarized 
proton beam, as a function of pion transverse momentum It- Theory curves are evaluated at 
XF = 0.4, = 200 GeV and A = 0.080 GeV. 
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TABLES 



TABLE L Partonic hard parts and corresponding color factors for quark-gluon and anti- 
quark-gluon (indicated by barred letters) subprocesses. Feynman diagrams are shown in Fig. 15. 
Cp, Cg, and Cg are color factors for spin-averaged, spin-dependent with initial-state interaction, 
and spin-dependent with final-state interaction subprocess, respectively. The explicit factors (-1) 
are due to the sign difference between quark and antiquark propagators with the same momentum. 
Calculations were done in Feynman gauge. 
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TABLE II. Partonic hard parts and corresponding color factors for subprocesses involving 



quarks and/or antiquarks. Feynman diagrams are shown in Fig. 17. In diagrams (i) and (j) both 
fermion arrows have been reversed, relative to (i) and (j). Cq, C^, and are color factors for 
spin-averaged, spin-dependent with initial-state interaction, and spin-dependent with final-state in- 
teraction, respectively. Flavor indices a and b correspond to the flavor of the quark (or antiquark) 
from the polarized hadron and unpolarized hadron, respectively, and c is the flavor of fragment- 
ing quark. The explicit factors (-1) are due to the sign difference between quark and antiquark 
propagators with the same momentum. Calculations were done in Feynman gauge. 
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